Resonant particle production with non-minimally coupled scalar fields 
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We investigate a resonant particle production of a scalar field \ coupled non-minimally to a 
spacetime curvature R (£Rx 2 ) as we ll as to an inflaton field <j) (g 2 4> 2 X 2 )- I n the case of g < 3 x 10~ 4 , 
£ effect assists p-resonance in certain parameter regimes. However, for g > 3 x 10~ 4 , g-resonance is 
not enhanced by £ effect because of £ suppression e ffect as well as a back reaction effect. If £ « — 4, 
the maximal fluctuation of produced x-particle is yj (x 2 ) max ~ 2 x f0 17 GeV for g < 1 x 10 -5 , which 
is larger than the minimally coupled case with g ~ f x 10~ 3 . 
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I. INTRODUCTION 



(N 
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In an inflationary scenario, most of the elementary particles in the Universe were created during the stage of 



reheating after inflation jlj . During the inflationary stage, an inflaton field slowly rolls down to a minimum of its 
potential. A reheating process turns on when the inflaton field begins to oscillate around the minimum of its potential. 
The original version of reheating scenario was first considered in the context of a new inflationary model in . The 



Oh! 
i 

particles are created through an interaction term between the inflaton and some fields. A phenomenological decay 



term to describe the fact that the inflaton field decays to other lighter particles (radiation) is included in the equation 



of the inflaton field 0, and the energy of the inflaton is transferred to their thermal energy. According to this scenario, 

_ _ i 

reheating temperature is determined by the decay rate T but not the initial value of <f>. The value of T is constrained 
to be small by the perturbation theory such as T < 10~ 20 AfpL, where Mpl is the Planck mass, so the reheating 
temperature is estimated as T r < 10 9 GeV. This temperature is not sufficient in order to produce baryon asymmetry 
based on Grand Unified Theories (GUTs). 

Recently, however, it has been recognized that reheating process begins by a parametric amplification of scalar par- 
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tides This initial evolutionary phase, which is called preheating stage, provides an explosive particle production 

and must be discussed separately from the perturbative decay of inflaton. There are many works about the preheating 
stage based on analytical investigations PHll|] as well as on numerical studies . The important feature with the 

existence of preheating stage is that the maximal value of produced fluctuation can be so large that it would result in 
a non-thermal phase transition p| and make baryogenesis at the GUT scale possible [jl6|, although the baryogenesis 
might be important in much lower energy scale, i.e. the electro- weak scale p7j . 

So far, we know two possible preheating scenarios: one is that the inflaton held 4> itself is transformed into many 
</>-particles from a coherent inflationary phase through a self- interaction such as X(j> , and the other is that another 
field, e.g. a scalar field \ is created through some coupling with the inflaton such as g 2 <t> 2 x 2 ■ Both cases were first 
discussed in S using Hartree mean-field approximation. For the former case, the preheating was studied ||| by 
making use of closed time path formalism (T^-|2^], since the preheating is an essentially non-equilibrium state. They 
analyzed the non-perturbative evolution of the inflaton fluctuations for the self-interacting massive inflaton by the 
method of the <3(7V)-vector model in the large N limit as well as by the Hartree factorization model, both of which are 
mean-field approximations. The 0(N)-vectoi model has an advantage to deal with the continuous symmetry while 
the Hartree factorization model is suitable to treat the discrete symmetry. To confirm the mean-field approximation, 
the fully non-linear numerical simulation including the scattering effect of created particles was also performed in 
|L2| for the simple A</> 4 model, finding that the variance (S(/) 2 ) ~ 10 _7 Mp L can be produced by the non-perturbative 
process of inflaton decay. They treated the scalar field fluctuations as classical ones, which would be justified because 
the fluctuations with rather low momenta are mainly produced by the parametric resonance. They found almost the 
same result as that by the mean-field approximation. 

As for the case only with a massive-inflaton potential (^m 2 (j> 2 ), the production of inflaton fluctuation will not be 
expected by a parametric resonance. Hence one usually introduces another scalar field x coupled to the inflaton field 
(f> through an interaction such as g 2 (j) 2 X 2 - For the production of 0-particles in self-interacting case, the coupling is 
constrained to be A ~ 10 -12 from the observation of COBE, while the coupling constant g in the present case is a 
free parameter. Hence we may find a higher production of ^-particles depending on the value of g. As for analytical 
investigation of x-field evolution, Kofman, Linde and Starobinsky [||J^] developed the consistent theory of preheating 
based on Mathieu equation for x-field. There are several numerical works devoted to the evolution of x-field by fully 
non-linear simulations p^JT^ |. A parametric resonance turns on from the broad resonance regime in certain values 



of the coupling g. The structure of resonance will be modified a lot in the expanding Universe compared with that 
in Minkowski space. The amplitude of the coherent (j> field decreases adiabatically because of the expansion of the 
Universe, and eventually the resonance terminates. In the case that the coupling g is small and the resonance band is 
narrow from the first stage of preheating, an efficient resonance will not be expected. For the large coupling constant, 
the resonance band is broad in the beginning, and then we find a considerable production of x-particles. In this case, 
X-field crosses many instability and stability bands even within one oscillation of the inflaton field, and stochastic 
resonance occurs. A significant amplification occurs in the broad resonance regime, but in certain values of coupling 
constant g, production via such a resonance will be suppressed by a back reaction of created particles. When such 
a back reaction is taken into account, a coherent oscillation of the inflaton field is broken due to the increase of the 
effective mass of inflaton, and the energy of %-field can not exceed much larger than that of i^-field at the final stage 
of preheating. Moreover, rescattering effect, which will also restrict x-particles production, becomes important at the 
final stage of preheating. 

Recently, another interesting mechanism for preheating called a geometric reheating model has been proposed by 
Bassett and Liberati pj| . They investigated the case that \ field is non-minimally coupled to spacetime curvature R 
for massive inflaton model and found that an explosive amplification of \ will be possible if the coupling constant £ is 
negative p3|. This is just because of unstable modes with a negative coupling [ p5[ . In this sense, it is different from 
ordinary parametric resonance. In particular, they studied how the homogeneous \ field is amplified and found that 
GUT scale gauge boson with mass ma ~ 10 16 GeV can be produced in preheating stage when £ is large negative. 

A natural question may arise: If we combine both effects of ^-resonance and ^-resonance, whether we find more 
effective production of %-particles, since small coupling g does not provide a sufficient production. It is of interest 
how the ^-resonance picture is modified by taking into account the effect of geometric reheating. If the assisting 
mechanism works and final abundance of \ becomes rather larger, it may change the ordinary preheating picture and 
affect the non-thermal phase transition and baryogenesis. We will discuss such a combined resonance in detail. We 
also study the evolution of \ field fluctuation summed by all possible momentum modes in detail including the back 
reaction effect, because the previous study is mainly devoted to the growth of zero momentum mode at the initial 
stage of preheating p3| . 

This paper is organized as follows. In the next section, we introduce basic equations in the non-minimal preheating 
that includes ordinary g-resonance. The results of g = 0, £ > case is presented in Sec. III. We study the analytical 
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estimation as well as numerical results, and compare them with each other. In Sec. IV, we analyze g = 0,£ < case. 
The structure of negative coupling instability is investigated. In Sec. V, the combined effect of g and £ resonance is 
studied. We present that at which values of g and £ parametric resonance becomes most efficient. Finally, we give 
our discussions and conclusions in the final section. 

II. BASIC EQUATIONS 



We consider a model that an inflaton field 4> is interacting with a scalar field %, which is non-minimally coupled 
with the spacetime curvature R, 



(2.1) 



where n 2 /8n = G = M p £ is Newton's gravitational constant, g and £ are coupling constants, and rn x is a mass of the 
X field. V((f>) is a potential of the inflaton field. In this paper we adopt the quadratic potential 

1 



v{4>) = 



(2.2) 



where m is a mass of the inflaton field, and we use the value of m — 1CP 6 A/pl that is obtained by fitting density 
perturbations to COBE data. 

Since non-minimal coupling between the spacetime curvature R and \ field makes the basic equations complicated, 



it is convenient to transform action (2.1) into the Einstein frame by a conformal transformation |2q |. We make the 
conformal transformation as follows: 



— O 2 



where f2 2 = 1 — £k 2 % 2 . The Lagrangian density in the Einstein frame becomes 



C 



where variables with a caret denote those in the Einstein frame, and 

2 _ 1 - (1 - 6flg«V 

To make the kinetic term of the x field canonical, we define a new scalar field X as 

X = [ F( X )d X , 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



by which the Lagrangian density is now 
C-- 



(2.7) 



Since we are interested in a preheating after inflation, as usual, we shall assume that the spacetime and the inflaton 
tj> give a classical background and the scalar field x is treated as a quantum field on that background. In the present 
model, however, there is some problem. We have performed a conformal transformation, which conformal factor fl 2 
includes quantum variable \ 2 ■ Then, in order not to discuss quantum gravity, the conformal factor should be replaced 
with some expectation value. Here, we regard the conformal factor f2 2 as 1 — r\ WAi where 



v^^ 2 (x 2 )- 



(2.8) 



(x 2 ) corresponds to the number density of the created x-particle, where (• • •) denotes an expectation value of some 
functional of \ ( or ^0- 

With such a transformation, we assume that the spacetime and the inflaton field <p are spatially homogeneous, and 
adopt the flat Friedmann-Robertson- Walker metric as the background spacetime; 



ds 2 



dt z + a 2 (t)dx 2 



(2.9) 



Hereafter, as we argue only in the Einstein frame, we drop a caret. The evolution of the scale factor a yields 









)'-t[ 


(a, 





1 



1 



mV + i<(VX) 2 ) + 



1 



_2(1-V) 2(1 — 77)2 T 2 U ;/ 2(1 - j]) 2 v x 
where a dot denotes a derivative with respect to time coordinate t. 
The evolution of the inflaton field is described by 



(2.10) 



V 



+ J—( m 2 +g 2 ( x 2 ))0 = Q. 
1 — 7] 



(2.11) 



a \ — r\ / 

Note that the fluctuation of the inflaton field 6(j) is not considered since we assume that the inflaton field is spatially 
homogeneous. However, the growth of the <f> fluctuation would be expected to appear as x field is amplified. Our 
investigations are limited in the sense that rescattering between x and S(f> fluctuation is not included. 



In Eq. (2.11), if we introduce a new scalar field (p such as 

where 6 3 = a 3 /(l — ry), then the field ip obeys the following equation: 
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(2.12) 



3/26 b 2 



1 — T) 



<p = 0. 



(2.13) 



Note that in the case of minimal coupling £ = 0, r\ vanishes, hence the coherent oscillation of ip is broken only by 



g 2 (x 2 ) term. As (x 2 ) grows, the effective mass of the infiaton 



"eff 



- g 2 (x 2 } gets large, i.e., oscillation becomes 



rapid. This effect, which is called back reaction effect to the infiaton field, suppresses the resonant particle creation. 
On the contrary, when £ effect is taken into account, 1 — rj is decreasing as (x 2 ) grows, hence one may expect this 
effect also changes the coherent oscillation of ip. However, as we will see later, this is not the case and we can neglect 
£ effect on the infiaton field in most cases. 



Next, let us consider the equation of X field. The Heisenberg equation of motion is derived from (2.7) 

0, 



x + 3-x - did'x - A 

a oX 



(2.14) 



where an index with a roman character denotes space coordinates. In order to study a quantum particle creation of 
X fields, we make the following mean field approximation with respect to X , which provides us a linearized equation 
for a quantum field X : 



X + 3-X - d l d l X - V 
a 



X = 0, 



(2.15) 



where V x = 2 V /(X 2 ) d/d{X 2 ). From the relation (§J), we have also assumed 



where (F 2 ) = [1 — (1 — 6£)7j]/(l — r\) 2 . Expanding the scalar fields X as 



(2.16) 



X = 



(27T) 



^ J (a k X k (t)e-* k -* + alX* k (ty k - x ) d 3 k, 



(2.17) 



where a k and a\ are the annihilation and creation operators, we find that X k obeys the following equation of motion: 



x k + s-x k + 

a 



- 2 +G{{x 2 )) 



X k = 0, 



(2.18) 



with 



G((X 2 )) 



(1 - t?) [1 - (1 - 6C)V] 

77 4+(l- 577)(l-6£) 



1-t] 1 - (1 - 6£)r/ 
The expectation values of X 2 and x 2 are expressed as 



(1 + 3r7)(m 2 . + g 2 <t> 2 ) + 2£« 2 m 2 (/> 2 - (1 - Z^n 2 ^ 2 

{(1 + f])(m 2 x + g 2 <t> 2 ) + 2^K 2 m 2 (f> 2 - (1 - ?y)^ 2 2 } 



(2.19) 



(X 2 ) = 2^2/ k 2 \X k \ 2 dk, (x 2 ) = 2^2/ k 2 \ Xk \ 2 dk. (2.20) 
Introducing the function Y k — a 3 ^ 2 X k , instead of X k , we find 

Y k + uj 2 Y k = 0, (2.21) 

where 

^ = K + G((x 2 ))-j(- + ^), (2.22) 
a z 4 \ a ar / 

which is a time dependent frequency of Y k . 

At the first stage of preheating when x _ P ar ticles are produced by quantum fluctuation, we find the occupation 

number in Y^-state by the Bogoliubov transformation as 

— t(wI"^)-5- <" 3 » 

However, after many ^-particles are created and each mode is amplified, the Xk field could be regarded as the classical 
field in a good approximation Jl3|,[l4|]. We may not have to use Eq. (2.23) at the classical stage. Khlebnikov and 



Tkachev developed the semiclassical description of fluctuation produced by the inflaton decay (T^j. In fact, they 
studied the inflaton decay by the classical equation of motion and performed a fully nonlinear calculation. To give 
the initial conditions, we follow their approach, i.e. the initial distribution for Y k is described as 



P[Y k ;t= 0] = ^L^exp [~2cu k (0)\Y k \ 2 ] , (2.24) 

and Yfc is correlated to Y k as 

Y k = -iw fc (0)Yfe. (2.25) 

We investigate the (x 2 ) evolution with those initial conditions as the semiclassical problem. 

From Eq. ( 2.19| ), we easily find that properties of the preheating in the non-minimal coupling theory are quite 



different from ordinary g-resonance. First, ^K 2 m 2 (j) 2 term and ^k 2 4> 2 term cause the resonance as well as interaction 
term g 2 <j) 2 with inflaton field. These different types of resonant terms either strengthen or weaken the resonance each 
other, depending on the coupling constants. Secondly, as the x-particles are produced significantly, the suppression 



effect by the second term in Eq. (2.19) becomes crucial. This means that in case of non-minimal coupling, we have 
to consider not only ordinary back reaction effects to the inflaton field and metric but also the suppression effect by 
this term. In what follows, we will investigate these issues in detail. 
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III. THE RESONANCE BY POSITIVE COUPLING £ 



In this section, we investigate g — 0, £ > case. First, however, we briefly review the ordinary g-resonance, i.e., 
the case with g ^ 0, £ = for comparison In this and next two sections, we mainly study the massless x held. 

For the massive case, we will give some discussion in the end of each section. 

If the back reaction of the \ ne ld to the inflaton field and metric are neglected, the inflaton field oscillates almost 
coherently with damping factor 3H and is approximately described as 



Q(t) sin mt. 



(3.1) 



The amplitude $(f) decreases with time as 



*(t) 



Mpi 



Znmt 



Then the time-dependent frequency of the each component Yk becomes 



2 "> i 2i2 ™ i 2/r2 • 2 j. 

uj h . = + g (f> = —t + q $ Sin mt. 



We can reduce Eq. (2.21) to the well known Mathieu equation 



where z = mt and 



dz 2 



[Ak-2qco82z]Y k = 0, 



(3.2) 



(3.3) 



(3.4) 



A k =2q 



(3.5) 



Am 2 



(3.6) 



k is normalized by m as k = fc/m. Stability or instability with Eq. (3.4) depends on the variables of Ak and q, 
which is shown by a stability-instability chart (see Fig. 1) (^Sj. In the unstable region (the lined region in Fig. 1), Yk 
grows exponentially as Yk oc exp(/ifcz) with the Floquet index \ik and particles with momentum k are produced. For 
small q, the width of the instability band is small and few fc-modes grow by this resonance. This is called the narrow 
resonance. On the other hand, for the large q, the resonance can occur for a broad range of the momentum fc-space. 
Since the growth rate fj, of the x _ P ar ticle becomes larger as the increase of the variable q, this resonance gives more 
efficient particle production than narrow one. This is called the broad resonance. Note that initial amplitude $ of 



the inflaton field and coupling constant g play important roles to decide whether resonance is narrow or broad. In 
the g-resonance case, the allowed region on the Mathieu chart is determined by the Eq (3.5) as 

A k > 2q. (3.7) 

Hence the broadest resonance is given by the limit line Ak — 2q. 

When q is sufficiently large initially, the resonance of each mode occurs stochastically [^|J^] . In this case, the frequency 
ujk decreases by cosmic expansion and Uk dramatically changes within each oscillation of the inflaton field, so the 
phases of \ field at successive moment of tj> = are not correlated each other. At the first stage of the resonance, the 
fields cross large number of instability bands. The periods when they are in the instability band are so short that 
the resonance can not occur efficiently compared with that in the Minkowski spacetime. However, nevertheless, the 
number of ^-particles can still grow exponentially. As q becomes smaller, the universe expansion slows down, and the 
fields stay in each resonance band for a longer time. When q drops down to about 1, the first instability band 

l-q-\q 2 <A k <l + q-^q\ (3.8) 



becomes important. When the variables decrease below the lower boundary of Eq. ( p.8| ) by the expansion of the 
Universe, the resonance terminates. We have to note here that there is another mechanism which terminates the 
resonance. When the initial value of q is large (g > 3.0 x 10~ 4 ), x-particles are produced efficiently and the back 
reaction onto the inflaton field cannot be ignored. This makes the oscillation of the inflaton field incoherent and 
finally stops the resonance. That is called the back reaction effect of the \ field. 

Now we turn to the case of g = 0, £ > 0. In this case the resonance occurs by the coupling between the x field and 



the spacetime curvature R |23|. The G((x 2 )) term defined by Eq. (2.19) becomes 

1 



(1 - rj) [1 - (1 - 60r?] 

77 4 +(1-60(1-5^) 
1-7? 1 - (1 - 6£)77 



2£K 2 m 2 2 - (1 - 37 ? )£k 2 (/> 2 



|2£k 2 to 2 2 - (1-77)$k 2 2 } 



(3.9) 



We easily find that if 77 increases up to the order of unity, G((% 2 )) diverges and the frequency of (2.22) increases to 
infinity. The resonance seems to continue effectively at the final stage of preheating. However, this is not the case. 
Our numerical calculation shows that the resonance terminates when 77 is much smaller than unity in any cases. With 
the condition of 77 <C 1, G((x 2 )) is rewritten as 



Note that there exists the suppression factor (1 + 6£t7) -2 , which makes the amplitude of the x field small effectively. 
Even in the case of r/ <C 1, we can expect that this suppression plays an important role when £ takes large values. 
We will see later that this effect becomes significant in the case of £ > 70. When r\ is small compared with unity, the 



back reaction to the inflaton field and the metric is negligible and we find from Eq. (2.13) that the (p field oscillates 
almost coherently. 

By using these relations, we can rewrite the equation of Y k in the form of the Mathieu equation 

d 2 Y k 



dz 2 

where 



[A k -2qcoa(2z-a)]Y k = 0, (3.11) 



* = 5 + 0OT' (3 - 12) 



y/(2 P ) 2 + 2(g K 2 $ 2 ) 2 

q= 2(1 + W 2 ' ( } 



Of 

4 
-1 



(V - $' 2 ) , (3.14) 



£k 2 $$' 



\2p + £k 2 ® 2 J v ' 

A prime denotes a derivative with respect to z. Eq. fl3.ll ) is not exactly the Mathieu equation, because it contains 



the time dependent phase term a. This term is due to the existence of (p 2 term in Eq. (3.10). If we define a new 
dimensionless time parameter t = mt/2-K, a is rewritten as 

1 ( -Ant \ 

" = tair ( l2n*P-l > (3 ' 16) 



where we used the relation (3.2). Since t represents the number of the oscillation of the inflaton field naively, a 
approaches zero rapidly within the first oscillation of the inflaton field. As a result, the time dependence of a does 
not affect the resonance process. 

Although the equation of Y k has the same form as the ordinary (7-resonance except for the time-dependent phase 
a, properties of resonance are not the same. After the first oscillation of the inflaton field, setting an initial time as 
t = 1/4 as Kofman et al. [||, we have the condition 

$>-. (3.17) 

m 
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Then the evolution of A/~ and q is approximated as: 



2 I 2 



4(1 + 6^) 



1 



(1 + 6£?7) 2 2tt 2 P' 



(3.18) 
(3.19) 



The line described by Eq. (3.18) lies below the one which is obtained in the ordinary g-resonance. As is studied in 
p3[ , the Ricci scalar can be replaced with (j) 2 term and non-minimal coupling provides another contribution. This is 



because the above relation (3.18) becomes different from the relation (3.5) of ordinary ^-resonance. In fact, we can 
easily estimate that without 4> 2 term, £ must be larger than 10 4 for the effective resonance, because we have only the 
similar coupling term to the (/-resonance in that case. Since the width of the instability bands are thick for large q (see 
Fig. 1), ^-resonance gives broader resonance and we may expect the efficient x-particle production. However, as is 



found from Eq. fl3.19|) , there are two factors which decrease the variable q. One is the decrease of <£> by the expansion 
of the Universe and the other is the suppression effect caused by the factor (1 + 6£,rj) ~ 2 , which appears only in the 
^-resonance. Which factor is more important depends on the parameter £ and the efficiency of the resonance, i.e., rj. 
When £i] < 0.1, this suppression effect can be neglected. Note that the back reaction effect is negligible because r\ is 
always small in any parameter. Hence, it may be worth stressing that the ^-resonance in the g = case will terminate 
only by passing through the resonance band, because the coherence of the inflaton is not broken. 

By the final value of q (=<?/) when the resonance ends, we can estimate the time tf when the resonance stops and 
the total amount of created x-partides (x 2 )/ = (x 2 )// m2 - I n the case of q > 1, particle creation occurs when the 
frequency u>k changes non-adiabatically, which condition is written as 



~dt' 



(3.20) 



Neglecting the <E> in calculation in the £k 2 (/> 2 term in Eq. ( 3.1C| ), cu 2 yields 

9 k 2 £k 2 to 2 



a 2 (l + 6^) ; 



(30 2 - $ 2 ) 



(3.21) 



and then the non-adiabatic condition (3.2C) becomes 



< 



3£« 2 m 2 



6£ 2 f]K 2 m 2 



(l + 6£7 ? ) 2 ^ (l + 6£?7) 3 



(3(f> 2 - $ 2 ) 



2/3 



^K 2 m 2 



(1 + e^y 



(3</> 2 - $ 2 



(3.22) 



Let us investigate the value of <j) when the r.h.s. in Eq. ( 3.22 ) takes the maximum value. Non-adiabatic amplification 
occurs mostly when <f> is passing through around the minimum of its potential (|(^>| <C Mpl), so we can set <fi ~ to$ in 
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Eq. 03.22D . Moreover, fj term is negligible compared with the former term, which is confirmed by numerical calculation. 
Then, Eq. ( [Oil can be approximately rewritten as 



< 



3£K 2 m 3 1 2/3 



{;K 2 m 2 



(1 + 6fr)S 



(30 2 - $ 2 ) . 



(3.23) 



Differentiating the r.h.s. of Eq (3.23) with respect to 0, we find that it takes the maximal value at 



1 / m$ 1 / m$ 

Was = Tpjl Y ~qTJ2 ~ 2 V Tj^J 2 ' 



(3.24) 



where C = 3^K 2 m 2 /(l + 6^) 2 . Then the maximum value of momentum yields 



C$ 2 
3m 2 



+ (2 4 / 3 - l) 



4m 



(3.25) 



where A:^^ = fc 2 aa3 ./m 2 . As we have the relation vC$/m = by Eq. ( 3.19 ), the maximum momentum is rewritten 
in terms of q as 



(3.26) 



This equation gives the maximum momentum for the resonance. While the resonance terminates when the variables 
Ak and q pass the curve of 



A k ^l-q-^q 2 , 



(3.27) 



in the Mathieu chart 28 1. Combining the relations ( 3.18 ), ( 3.26 ) and ( 3.27 ), we obtain the equation with respect to 
q/ as 



q 2 +24q f + 6^qJ+ 8 = 0, 

resulting qf = 0.2165 ~ 1/5. If we adopt the typical momentum fc*, which would be defined as 

h 

, ""max 

in stead of k max , we find 

3^ + 56^ + 9^/97 + 24 = 0, 



(3.28) 



(3.29) 



(3.30) 



and q f = 0.3353 » 1/3. 
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The analytic expression of Eq. ( 3.28D or ( |3.30 ) is only an approximation because the maximal momentum k n 



will be changed when q drops down to q < 1. However, it gives good agreement with numerical calculation. First, 
we examine the evolution of (x 2 )j comparing the results of numerical calculation with the analytical estimation. For 
£ < 10, numerical calculation shows no effective production of (x 2 )- This means that the initial value of q is so small 
and it decreases so fast because of the expansion of the Universe that the x - field does not stay in the instability bands 
for enough time. For 10 < £ < 70, the resonance occurs because the period during which the x field stays in the broad 
resonance band becomes longer. For example, in the case of £ = 50, the value of (x 2 ) first increases exponentially by 
the passage of time and reaches to its maximum value (x 2 ) / = 3.090 x 10 3 at it = 2.85 (Fig. 2(a)). Although (x 2 ) / 
is not large enough for the efficient preheating, the parametric resonance evidently occurs. Since the maximal value 



of £77 is £77/ ps 1.9 x 10~ 4 , we can ignore the suppression factor (1 + 6(,r/)~ 2 in Eq. ( 3.19 ). Hence the resonance is 
terminated by the expansion of the Universe. The final value of q obtained by the numerical calculation is qf ~ 0.312, 
which is almost the same as the analytically estimated value with the typical momentum fc*. In the case that the 
suppression effect is neglected, we can estimate the time when the resonance ceases by 



For the case of £ = 50, it = 2.76, which is close to the numerical value it = 2.85. After (x 2 ) reaches its maximum 
value, it decreases monotonically. This is the adiabatic damping due to the expansion of the Universe. 

In order to see the ^-dependence of the numbers of created particles, we depict (x 2 ) / in terms of £ in Fig. 3. It shows 
that although the created particle first increases as the coupling constant £ gets large, it rather decreases beyond a 
critical value £ c , which can be understood as follows. 

When £ is less than about 100, (x 2 ) / increases as the coupling constant £ gets larger. This is just because the 
initial value of q is larger and then the resonance begins in the broader bands. For 70 < £ < 200, the suppression 
factor becomes important. The £ = 100 case is shown in Fig. 2(b). At the first stage, (x 2 ) increases rapidly with 
the larger growth rate than the case with £ = 50, but it reaches its maximal value (x 2 )/ = 9.550 x 10 5 soon. Since 
£77/ = 0.240 at the maximum point, the suppression effect by (1 + 6^77) ~ 2 cannot be ignored. In Fig. 3, we find that 
(x 2 )/ is almost flat around ~ 10 6 in the parameter range £ = 100 ~ 200. 

For £ > 200, the suppression effect by large £ is more effective than increasing of (x 2 )- At the first stage of 
preheating, the growth rate of (x 2 ) becomes larger as the increase of £ (Fig. 2(c)). However, the resonance soon 
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terminates with the smaller final value (x 2 ) / than in the case of £ « 100 by the suppression effect. For example, in 
the case with £ = 1000, (x 2 ) / = 8.318 x 10 4 and £77/ = 2.090 at if = 1.31, which is much smaller that the amount in 
the case with £ = 100. In fact, (x 2 ) / decreases monotonically by the suppression effect beyond £ c ~ 100-200 (Fig. 3). 
In the case of £ > 200, we can estimate the value of (x 2 )/ as follows. Assuming that the mode with the typical 
momentum fc* is the leading mode of the growth of (x 2 )> we can rewrite Eq. ( 3.1 9| ) as 



{x)f ~ 48^ 2 2^Wf ~ ) ' (3 ' 3 } 

at the maximal point. In our numerical calculation, we find that if is well approximated by the constant value 1.31 
when £ is greater than 500. Then, (x 2 )/ decreases as 

(f)/«f 3/2 , (3.33) 

which is confirmed by our numerical result given in Fig. 3 for £ > 500. 

In Table I, we also show our numerical results and the estimated values ( 3.32] ) for the maximal value (x 2 ) /■ The 



present analytical estimation gives a good approximation to the numerical results. A small difference between those 



may be due to the naive condition of Eq. (3.20). If we take into account all momenta larger than fc*, we will obtain 
the larger values of (x 2 ) analytic-, which estimation gives closer value to the numerical one. 

Finally, we should mention the effect of a mass of the x-held. In general, the mass effect of the x-field works as a 



suppression factor, because the relation between Ak and q in (3.12) is modified as 



k 2 2n ml 

a z (1 + be;?/) m 

However, since the resonance band is wider in ^-resonance compared with (/-resonance, the effective production of 
X-particles is still expected even if the mass term is taken into account. In fact we have found that if the mass is order 
10 13 GeV (namely the same order of the inflaton mass), x-particles are created as same as the massless case when £ is 
about 100 ~ 200. The example is given in Fig. 4. The final value of (x 2 ) in the case of g = 0, £ = 200, m x = m is about 
10 6 , which is almost the same of the massless case. However, for £ < 1000, it is difficult to have a resonant production 
of x-particles with the mass larger than 100m(~ 10 15 GeV), because the relation between Ak and q deviates from 
the resonance bands due to the mass term. When £ is much larger than 1000, although it is not likely, we still expect 
production of massive x _ P ar ticle whose mass is of order 10 15 GeV at the initial stage of preheating. However, as the 
production of x-particle proceeds, £77 suppression becomes efficient for such a large value of £, resulting in a small 
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total amount of (x 2 ). Namely, massive %-particle production is possible when £ is sufficiently large, but the final 
amount of x-particle gets small by the £77 suppression effect. Whether or not such a small amount of production can 
still provide us the baryonsynthesis is another problem, which we do not investigate here. 

IV. THE RESONANCE BY NEGATIVE COUPLING £ 

Next, we investigate the case of g = 0, £ < 0. Since 77 is small compared with unity as we will see later by numerical 



calculation, we can use the Mathieu equation (3.11) ~ (3.15) to analyze our numerical results in this case as well 



Neglecting $ term, the relation between Ak and q is now 

2 P 

A k n-~q + —, (4.1) 

O (X 

where 

qK I £L (42) 

Note that Ak can take a negative value because of the negative coupling £. This fact makes the properties of resonance 
quite different from the case of £ > as was pointed out in Ref. f23|| . In the regions of Ak < 0, a new instability band 
(zeroth instability band) extends below some curve which is approximated by Ak ~ —q 2 /2 when q is small (see Fig. 1). 
One of the important features is that this zeroth instability band reaches Ak — in the limit of q = 0. Consider the 
modes with very small momentum k. The line of Ak = — 2g/3 crosses to the curve Ak » — g 2 /2 at q ss 1.4. Hence, 
even if q evolves below unity by the expansion of the Universe, there remain some unstable modes in this instability 
band. Since the Floquet index jik of this instability band behaves as fj,k ~ ^fq when q is small, resonance continues 
to occur until q vanishes. Moreover, most of the Ak < region is covered by either the zeroth or higher instability 
bands. Almost all modes contribute to the resonant process in almost all the time with high growth rate. Hence one 
can expect the very effective particle production. 

We show numerical results for the evolution of (\ 2 ) in Fig. 5. For £ = —20, q takes initially the value ~ 16 and 
modes are either in the first or in higher instability bands, depending on the momentum. However, after the first 
oscillation of the inflaton field, q decreases to ~ 1 and some modes enter into the zeroth instability band, (x 2 ) 
increases exponentially mainly by those modes during the first several oscillations of the inflaton field. Thereafter the 
growth rate becomes small and (x 2 ) approaches the constant value (x 2 ) / = 4.1 x 10 7 . At this stage, the production 
of x-particle balances with the dilution by the Hubble expansion. For £ = —100, (x 2 ) increases more rapidly than the 
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case of £ = —20 at the first stage and reaches to its maximum value (x 2 ) / = 4.4 x 10 6 at t = 7.7. The main growth 
of (x 2 ) occurs until i = 2, when q is of order 1. After i = 2, the x-field enters the zeroth instability band, and the 
increase of (x 2 ) stops at i = 7.7. After that, (x 2 ) takes almost a constant value. This behavior is universal for the 
case with £ < 0. To understand the present results more deeply, we shall estimate the value of (x 2 )f- 
From the relations fl2.16| ) and ( |2.20D , we can obtain the following relation: 

d 



dt 



(x 2 



1 (x 2 ) i 

(F 2 )(X 2 ) a 3 \dt 



:<X')- 3H(Y 2 ) , 



! (x 2 )(x 2 



(2(1 - 3H) 



(4.3) 



where we have used the expression (Y 2 ) = a 3 (X 2 ) ~ e 2 ^ mt with /i being the Floquet index and 5 = H/m. The 
growth of (x 2 ) stops when the r.h.s. vanishes, i.e., when the dilution effect by expansion of the Universe surpasses 
the particle creation rate. Neglecting the back reaction effect on the metric, the evolution of the Hubble parameter is 
approximately written by 



H 



I Aw $ 



f 



3 Mpl 3nt 

While, the relation between ji and A in the zeroth instability band can be written by |26 



(4.4) 



-A- 



(A-l)q 2 



1/2 



2{A- l) 2 -q 2 

Considering the modes which are close to the line of A = —2q/3, the Floquet index jji is given as 



(4.5) 



2 1 ~ 1 
^ ~ V T ~ 1 + 6£?7 V 3 7rt ' 



(4.6) 



Until (x 2 ) reaches to its maximum value, \i decreases faster than 1/i by the existence of the suppression factor 
1/(1 + 6£?/). Because of this behavior, the Hubble expansion, which decreases as ~ 1/i, will catch up with the 
inflaton decay, and then the preheating is terminated. Namely, the main factor which stops the growth of (x 2 ) is the 



production of (x ) itself. Substituting Eqs. (4.4), (4.6) to Eq. (4.3), the growth rate of (x ) is approximately given as 



dV X ' 



'(x 2 )(x 2 



(F 2 



TTt 



1 



1 + 6^ 



- 1 



(4.7) 



When the term in the square bracket vanishes, the growth of (x 2 ) ceases. By this condition we find the final value of 

(X 2 > as 

(M PL /m) 2 



(x 2 )/ 



48tt£ 2 



3 



(4.8) 



1G 



Note that the final value of (x 2 ) depends on only £, but not on if. Once (x 2 ) reaches to its maximum, n and H both 
decreases as 1/i. Note that in the case of g = 0, £ > 0, the dependence of fi is /i ~ q ~ 1/P when q is close to zero 
and that (x 2 ) does not approach to constant value but adiabatically decreases by the Hubble expansion. 



From Eq. (4.7), we find that the resonance does not occur in the case of < |£| < 0.75, and this is confirmed by 
numerical calculation. In this parameter range, the growth rate \i is small as compared with the Hubble expansion 



rate. By Eq. (4.8), we guess that (x 2 )/ may take the maximal value at £ « —1.33. However, this is not the case 
actually. In —3 < £ < — 1 case, since the creation rate of %-particle is small, it takes much time to complete the 
X-particle production. As time passes, a contribution of produced x-particles to the Hubble expansion rate becomes 



comparable to that of the inflaton field, and the estimation (4.8) can not be applied. The growth of (x 2 ) stops before 



£77 suppression effect in Eq. (4.7) becomes significant. Hence for —3 < £ < — 1, although the x _ P a rticle production 
is possible, the final abundance of ^-particle is not so large. For example, in £ = —2 case, the numerical value of 
the final abundance is (x 2 )/ = 3.2 x 10 3 , which is much smaller than the estimated value (x 2 )/ = 1-0 x 10 9 by 
Eq. (|4.Sj). On the other hand, for £ < — 3, £77 suppression effect plays a crucial role to teriminate the increase of (x 2 ). 
Numerically, (x 2 )/ takes maximal value (x 2 )max ~ 3.2 x 10 8 at if ~ 3.3 x 10 5 when £ ~ —4. In this case, although 
the growth rate is still small compared with other cases as g > 0, £ = and g — 0, £ > 0, the final fluctuation is large. 
For £ < — 4, (X 2 ) continues to grow until x-particles are significantly produced, and finally £77 suppression effect 



terminates the resonance. In these cases, the analytic estimation based on Eq. (4.8) can be applied. By Eq. (4.8), we 
expect that (x 2 ) / decreases as the decrease of £ (< — 4). For example, numerical values of the final fluctuation are 
(X 2 )/ = 1.122 x 10 8 for £ = -10 and (x 2 )/ = 4.130 x 10 7 for £ = -20. In £ = -10 case, if is about if - 1.0 x 10 3 
and the growth rate is /i ~ 0.001. For £ = —20, the growth rate increases up to \i ~ 0.1, and the resonance ends at 
if = 28.28. In the case of |£| 3> 1, the resonant time becomes very short by this suppression and the final value of 



(x 2 ) is reduced. Eq. (4.8) shows that (x 2 )/ decreases as |£| -3 / 2 f° r \C\ ^ 1- 

In TABLE II, we show the analytically estimated and numerical values of (x 2 ) / for various cases. We find that the 
analytical estimation gives good agreement for £ < — 4. The small discrepancy comes from the fact that we mainly 
considered the modes close to k — 0, which gives the largest contribution in most of the stages. Also, since the actual 



Hubble parameter is larger than is estimated by Eq. (4.4), this decreases the estimated value of (x 2 ) /■ Note that the 
final value of (x 2 ) / = 3.162 x 10 8 at £ — —4 is larger than the maximal value in the case of g — 0, £ > ((x 2 )max ~ 10 6 
at £ = 100 ~ 200). Since the resonance bands of negative coupling are broader than those of positive coupling in 
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previous section, the more x-particles are created and the final value of (x 2 ) becomes larger. However, when the value 
of |£| becomes large, (x 2 )/ decreases as |£|~ 3 / 2 , which is the same as positive coupling case. We conclude that the 
suppression effect controls the finial value of (x 2 ) for large |£| in ^-resonance case. 

When the mass of x-held is taken into account, the resonance is suppressed as the case of g = 0,£ > 0. However, 
since the resonance band is broader than those of g = 0, £ > case, we expect considerable production of x-particles 
for £ < 0. One important property in negative £ case is that massive %-particle is hard to be created once x-field 
enters the zeroth instability band. This is because x-field deviates from the zeroth instability region by the mass effect. 
Even in the case of £ = —20 ~ —30 when massless x-particles are effectively produced as (x 2 )/ > 5 x 10 7 , massive 
X-particle production is strongly suppressed because x-held deviates from the zeroth instability band soon after the 
first oscillation of the inflaton field. For example, in the case of g = 0,£ = — 30,to x = m, we find (x 2 )/ ~ 10 3 , which 
is not an effective resonance. Rather, in the case of g = 0, £ = —100 ~ —200, in spite of the £ry suppression effect, the 
produced ^-particles exceed more than (x 2 ) / = 10 6 for m x — m (Fig. 6). However, if m x > 10m, it becomes difficult 
to produce more than (X 2 ) / ~ 10 6 even for £ = —100 ~ —200. We show in Fig. 7 the final abundance of x-particle 
as a function of m x for £ = —1000 case. In this case, we find that x-particle whose mass is more than m x ~ 100m 
can not be created. In order to produce the GUT scale particles (m x ~ 10 3 m ~ 10 16 GeV), the value of Aj- is more 
than Ak — 10 6 since m 2 /m 2 ~ 10 6 in Eq. ( |3.34| ). This means that the resonance does not happen at all from the 
very beginning unless £ < — 10 6 [Note q « |£| initially]. Moreover, when |£| is extremely large such as £ ~ — 10 6 , 
X-particle with the mass of order 10 16 GeV can be created initially, but the final amount will be largely reduced by 
the £?7 suppression effect. As is the same with g = 0, £ > case, whether or not such a small amount of production 
can provide us the baryonsynthesis is another problem. We then conclude that for massive case, effective x-particle 
production is expected when the \ ne ld does not deviate from instability bands initially and (j\ suppression effect is 
not too strong. 

V. THE COMBINED RESONANCE 

In previous two sections, we have shown that the non-minimal coupling will give rise to the parametric resonance 
depending on the value of coupling constant £. Here we shall investigate a combination of ^-resonance and ^-resonance 
(the combined resonance). When g ^ 0, g 2 (x 2 } term in the equation for the inflaton field ( 2.13] ) becomes important 



at the final stage of the evolution and changes the effective mass of the inflaton field. By this back reaction effect, 
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the production of ^-particles i s suppressed. Similarly, produced x-particles will affect the metric by the coupling 
between the x an d inflaton fields. However, first, in order to examine the naive structure of the present system, we 



shall rewrite the equation for the \ field under the condition of g 2 (x 2 ) "C m 2 and t)« 1. With Eqs. (3.1) and (|3.2|), 
the equation for the \ field is reduced to the form of Mathieu equation (3.11) with ( |3.12| ), ( 3.15| ), and q and p being 
defined as 



q = 



y/(2p) 2 + ( 5 2 $ 2 /m 2 + 2£k 2 $ 2 )£k 2 $ 2 
2(1 + 6^7 7 ) 2 



1 (9 * + 



(5.1) 
(5.2) 



With the condition (3.17), Ak and q are rewritten as 



fc 2 g 2 $ 2 /m 2 + £« 2 $ 2 _ k 2 5 2 + 8<(m/M PL ) 2 

fe ~ ^2 + 9| ff2$ 2/ TO 2 + 3 ^2 $ 2| w ^2 + « 



4(1 + 6^77) 2 



5 2 $ 2 



3£k 2 $ : 



-|5 2 + 247re(m/M PL ) 2 | 
(Af PL /m) 2 
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(5.3) 
(5.4) 



48tt 3 (1 + 6£r)) 2 P 

On the Mathieu chart, there are two quantities which determine the efficiency of the resonance. One is the initial 
value qi. This determines the Floquet index, i.e., the growth rate of the x-particle. It becomes large as g and £ get large 
as seen from Eq. ( |5.4|) . In general, if we can neglect the back reaction effect, we expect more efficient resonance for large 
qi . Another quantity is the gradient of the line on which the variables Ak and q trace during evolution of the system. 
It determines how many modes contribute to resonance. In the present case, the gradient is given by the ratio of £ to 
g 2 . In the minimal coupling case, £ = 0, so the typical line is Ak = 2q. In the non-minimal case, however, there are 
several types of lines which show different behaviors. Then we shall first classify the Mathieu chart into several regions 
by changing £/<? 2 from oo to — oo (see Fig. 8). When £/<? 2 = oo, the resonance is the similar to the positive ^-resonance 
discussed in Sec. III. When < £/<? 2 < oo, the resonance occurs below the Ak — 2q line. We call the resonance in 
this parameter region "new-broad resonance". When £ = 0, resonance occurs only by g-coupling: Ak = 2q + k 2 /a 2 
and the ordinary g-resonance is recovered. If £ becomes negative, Ak and q enter the region of Ak > 2q, and we call 
the resonance in this parameter region "ordinary-broad resonance". When £/<? 2 = — (AfpL/m) 2 /247r, q vanishes and 
no resonance occurs. As £,/g 2 decreases further, the gradient gets smaller and the resonance becomes broader again. 
When £/.g 2 = — (Mpl/to) 2 /167t, the typical line is Ak — 2q, which is same as the ordinary ^-resonance in spite of 
the existence of ^-coupling. The parameter range — (Mpl/to) 2 /127t < £/<? 2 < — (Mp^ / m) 2 / 16ir , corresponds to the 
region between Ak = 2q and Ak — 2g/3, and the system enters the new-broad resonance region again. As £/ 1 g 2 gets 
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smaller further, Ak and q enter the region of Ak < 2q/3. We call the resonance in this region "wide resonance". 
When £/g 2 < — (MpL/m) 2 /87r, the gradient becomes negative. £/g 2 = — oo corresponds to the negative ^-resonance 
discussed in the previous section. In what follows, we will examine each cases one by one in detail paying attention 
to which parameter range in ^-coupling assists the (/-resonance. 

Before proceeding the individual investigation, we shall give some criteria for the back reaction effect on the inflaton 
field and for the suppression effect by £?y term. 

1. For the back reaction effect on the inflaton field, we shall adopt the criterion 

9 2 (X 2 ) < 0.1. (5.5) 

If the (x 2 ) exceeds this criterion, the back reaction becomes important. In the case of g ^ 0, £ = 0, our numerical 
analysis shows that this criterion corresponds to the condition of g < 3 x 10~ 4 . 

2. For the £77 suppression effect, we adopt the criterion 

£77 < 0.1. (5.6) 

We will see that those effects are crucial for the effective x-particle creation at the final stage of preheating. 

A. £ > (new-broad resonance) 

This case corresponds the parameter range, 

2 fc 2 fc 2 

- q+ -<A k <2 q+ -,, (5.7) 

on the Mathieu chart (Fig. 8(a)) and the resonance band is broad compared with usual (/-resonance case. We shall 
discuss the following three cases separately. 

1. j«3x 1CT 4 case 

In the case of £ = 0, i.e. the ordinary (j-resonance, the width of the instability bands are small when j<3x 10~ 4 . 
The x-particles are actually created only in the first instability band, and the resonance occurs in the narrow band 
from the beginning. As a result, x-particle is not produced efficiently. By taking ^-coupling into account, however, 
the resonance is assisted since initial value of q becomes larger and the resonance band becomes broader as seen from 
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Eqs. ( |5.3| ) and (5.4). In order to obtain (x 2 )/ ~ 10 6 , we need the coupling constant £ = 100 ~ 200. Because the 
resonance occurs mainly by ^-coupling, the properties of resonance is almost the same as g = 0, £ > case. For fixed 
value £, we find that g makes the resonance band narrower and the resonance less effective, although it makes the 
initial value of q larger. In fact, the numerical calculation shows that (x 2 ) f with g is smaller than that without g and 
its growth rate of (x 2 ) is also smaller. In this sense, (/-coupling weakens the ^-resonance. In this parameter region, 
the back reaction effect on the inflaton mass is negligible and the suppression by £77 term terminates the x-particle 
production. 

2. j~3x 10~ 4 case 

In this case, we find from Fig. 9 that the resonance is divided into two stages as in the ordinary g-resonance, namely 
first broad resonance regime and following narrow resonance regime. We put the suffix "6" to the variables just after 
the broad resonance regime as ib and q^. Although the resonance starts from the broad resonance region even in the 
£ = case, the resonance band becomes broader by adding the ^-coupling. Furthermore, the initial value of q becomes 



large as seen in Eq. (5.4). Since more modes contribute to the resonance, both the growth rate of (x 2 ) and the value 
of (x 2 )t> get larger (see Fig. 9 and Table III). This means that the ^-coupling supports the (/-resonance well in the 
broad resonance regime. After the resonance enters the narrow resonance regime (t w 8), the growth rate hardly 
depends on £. At the final stage, however, we cannot neglect the suppression effect by the £77 term. As £ becomes 
large, the variables Ak and q cross the lower boundary of the first instability band at an earlier time and the final 
value of (x 2 ) / is suppressed a little compared with that in the ordinary g-resonance. In this case, the back reaction 
effect is marginally less important. 

3. g > 3 x 1CT 4 case 

In this parameter region, the back reaction to the inflaton field should be taken into consideration. In the case 
of £ = 0, x-particles are produced quite effectively in the broad resonance regimes and (x 2 )/ takes maximal value 
(x 2 )max ~ 5.0 x 10 7 for g sw 1 x 10~ 3 , in which case the initial value of q (^ w 1.075 x 10 4 ) is large enough as 



estimated by Eq. (5.4). As g increases, the back reaction term g 2 (x 2 ) gives significant effect even small value of 
(x 2 ) and (x 2 )/ rather decreases. Even if we take ^-coupling into account, it does not change the growth rate at the 
first stage of the resonance very much when g is large (see Fig. 10). For example, in the case of g — 1.0 x 10 3 and 
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£ = 100, the initial value estimated by Eq. (5.4) is q L w 1.083 x 10 4 . This is almost same as the case of £ = 0. 
Furthermore, the relation between Ak and q is Ak ~ 1.99g+ fc 2 /a 2 , which shows also the same broadness. Hence the 
initial stage of the resonance is governed by the g-coupling. In order to find the effect of the ^-coupling, we may set £ 
very large (£ > 10 4 for g = 1 X 10~ 3 ). Such a large £ coupling, however, causes extremely strong suppression and we 
do not expect the efficient resonance. In the final stage, either the suppression effect by £77 or the back reaction effect 
terminates the resonance depending on the coupling constants g and £. Generally speaking, when £ takes a large 
value, the suppression factor becomes important before the created ^-particles cause the significant back reaction. 



By two criteria (5.5) and (5J3), we can estimate a rough dependence. Once the created particles exceed one of the 
criteria first, the resonance will terminate. Since the r.h.s. of these criteria are the same, comparing the l.h.s. terms, 
two effects become equivalent when 



.9 2 = Ml^U 2 - (5-8) 



s M PL 

Hence, if g 2 > 87t(to/Mpl) 2 £ 2 , the back reaction terminates the resonance. Otherwise, the resonance is terminated 
by the suppression effect. This estimation has been confirmed by the numerical calculations. As a result, it is difficult 
to create %-particles more than (x 2 ) / ~ 5.0 x 10 7 even if ^-coupling is taken into account. 

B. — (Mph/m) 2 /16ty < ^/g 2 < (ordinary-broad resonance) 



As we can see by Eq. ( |5.3| ), both terms of g and £ suppress the instability in this case. The gradient of the q-Ak 
line becomes steep, and this parameter range corresponds to the region between the line of Ak = 2q + k 2 /a 2 and Ak 
axis in the Mathieu chart. This means that the resonance occurs in narrower bands than that in the case of £ = 0. 
In particular, for £,/g 2 — — (Mp^/m) 2 /24ir, the variable q vanishes. Since the width of any instability bands vanishes 
for Ak > 0, the resonance does not occur at all although both ^-coupling and ^-coupling exist. As for the initial value 



of q, both terms of g and £ suppress it as seen from Eq. (5.4). As a result ^-coupling always suppresses g-resonance 



in any parameter range. This effect is remarkable when q is small. However, for large g, £ terms is negligible because 
both the gradient of the q-Ak line and qi depend on square of g while they depend linearly on £. Hence the properties 
of the resonance are almost the same as those in ordinary g-resonance if g > 3.0 x 10~ 4 and < |£| < 100. 
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C. -(Mp L /m) 2 /127r < £/g 2 < -(Af PL /m) 2 /167r (new-broad resonance) 



This parameter range corresponds to the regions between Ak = 2q + k 2 /a 2 and Ak = 2q/3 + k 2 /a 2 on the Mathieu 
chart (Fig. 8(b)). The resonance gives the same broadness as the case discussed in the subsection A. However, there 
is a different point in an efficiency of the resonance. In the present case, q evolves in the range of 



(1 + 6^) 2 67r 2 t 2 * (1 + 6£?7) 2 47r 2 ^ 2 ' 



Comparing this with the Eq. (3.19) in the case of g = 0, £ > 0, q in Eq. (5.E) is smaller than that in the case of 
g = by factor 2-3 for the fixed value of |£|. Remembering that for g = 0,f > the value of £ must be £ > 100 
to achieve (x 2 )/ ~ 10 6 , we can estimate that |£| is needed at least |£| > 200 for the efficient resonance in the case 
°f C/s 2 ~ — (Mpl / m) 2 / 12tt . This means that the parameters which cause the effective resonance are a little more 
restrictive than the f -resonance discussed in Sec. III. In the case of f/<? 2 ~ — (Mp^ / in) 2 / 16ir , q evolves on the line 
A). w 2q + k 2 /a 2 and the broadness of the resonance bands is similar to the ordinary g-resonance. For this reason, 
we need the value of |£| that is more than 10000 for the effective resonance, but such a large |£| causes the strong £?y 
suppression and the final value of (x 2 ) is reduced. We show the numerical results on the Ak = 2g/3 line in Table IV, 
have confirmed the above analysis. If |£| is less than 100, the resonance hardly occurs. For example, when £ = —100 
and g — 6.140 x 10~ 5 , the final value of (x 2 ) is (x 2 )/ = 5-0 x 10. On the other hand, in the case of £ = —200 and 
g = 7.089 x 10~ 5 , the final value is (x 2 )/ = 1-0 x 10 5 . ^-resonance is in fact assisted by ^-coupling in this case. 
However, it is difficult to produce the x _ particles more than (x 2 )/ ~ 10 6 in this parameter range. 

D. —00 < £/g 2 < — (AfpL/m) 2 /127r (wide resonance) 

This parameter range corresponds to the regions between Ak — 2q/3 + k 2 /a 2 and Ak = — 2q/3 + k 2 /a 2 (Fig. 8(b)). 
Since the contribution from the negative ^-coupling surpasses the (/-coupling, the resonance band becomes quite broad 
as £/g 2 decreases. For — (Mp-^/m) 2 /%tt < £/g 2 < — (Mpl/to) 2 /127t, the gradient of the q-Ak line is positive and the 
resonance terminates when the variables q and Ak cross the lower boundary of the first instability band as the former 
cases. After the resonance terminates, (x 2 ) decreases gradually by the adiabatic expansion. We show the result of 
the numerical calculations in the £/g 2 = — (Mp^ / m) 2 / 8tt case, which correspond to Ak = 0, in Table V. In the case 
of |f I < 50, resonance does not occur because the initial value of q is too small. However as |£| becomes large, (x 2 )/ 
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increases and for £ = —200 and g — 7.089 x 10~ 5 , (x 2 )/ = 1-0 X 10 6 . From these results we can conclude that 
resonance is more effective than the previous case. 

When £/g 2 < — (Mpl / m) 2 / '8ir , the q-Ak line passes through the zeroth instability band near its boundary. Since the 
Floquet index is small there, it takes long time to reach the finial value of (x 2 ) even for the large |£| . (x 2 ) / is determined 



by the balance between the creation rate of the x-particle and the expansion rate of the Universe (Eq. (4.8)). As 
£/<? 2 approaches to — oo, the properties of the resonance are almost the same as the case of g = 0, £ < 0, and then 
^-coupling assists g-resonance. Since the existence of (/-coupling makes the growth rate small for fixed value of £, 
resonance is most efficient in the case of g = 0. 

E. summary 

We shall summarize the properties of the resonance and the suppression effect. On the £-g diagram (Fig. 11), we 
show which suppression factor is significant in various coupling regimes. We can easily find that the back reaction 
effect is dominant in the case of g > 3 x 10~ 4 and g > 5 x 10~ 6 |£| (the lined region in Fig. 11). However, in 
the case of g < 3 x 10~ 4 and |£| 3> 1, the £77 suppression effect becomes more important. In fact, £77 suppression 
appears either for £ > 100 or £/<? 2 < — 0.1(Mpl/to) 2 (the shaded region in Fig. 11). In the parameter range like 
— (AfpL/m) 2 /127r < £/g 2 < — (Mp^ / m) 2 / 16ir (new-broad resonance with £ < case), however, although |£| is of order 
100, the effective x-particle production will not be expected. In the region where the back reaction or suppression 
effect becomes significant, we may expect a large amount of particle production, although it will be reduced for 
extremely large |£|. 

In Fig. 12, we show (x 2 )/ in terms of £ and g in 3-dimensional plane. From this figure, we find three plateaus; one 
corresponds to the back reaction region with large g, and the other two correspond to £77 suppression regions shown 
in Fig. 11. (x 2 ) / takes maximal value of (x 2 )max ~ 3.2 x 10 8 at g < 1 X 10~ 5 , £ w —A. In other regions, we find that 
the numerical result agrees well with our analysis. 

In this section, we have studied the combined resonance of g and £. For g > 3 x 10~ 4 , when the back reaction 
effect is important, the final values of (x 2 ) does not become larger than those of ^-resonance significantly even if 
^-coupling is taken into account. This is due to both of the £?y suppression and the back reaction. In the case that 
|£| is large as |£| > 100, the £77 suppression effect becomes more important than the back reaction. For g < 3 x 10 -4 , 
(7-resonance is sometimes assisted by ^-coupling. For example, for 100 < £ < 200 and —00 < £/<? 2 < — (Mpl/to) 2 /167t, 
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^-coupling assists ^-resonance. In particular, for g < 1 x 10 5 , the resonance structure is almost the same as the case 
of g = 0, £ 7^ 0, which means that it is essentially the ^-resonance. 

VI. CONCLUSIONS AND DISCUSSIONS 

In this paper we have examined the properties of resonance with non-minimally coupled scalar field \ m preheating 
phase. We have found that effective resonance is possible only by a non-minimal coupling £i?x 2 in a certain range of 
parameter £. In the case of g — 0,£ > 0, the relation of Ak and q in the Mathieu chart is Ak = 2q/3 + k 2 /a 2 and 
the resonance band becomes broader in comparison with g-resonance (Ak = 2q + k 2 / a 2 ). This is due to the existence 
of £,K 2 (f> 2 term and the structure of resonance is different from that of g-resonance. Without £k 2 </> 2 term, £ must be 
larger than 10 4 for the effective resonance. However, as we have shown here, the effective resonance is possible for 
£ ~ 0(100). For example in the case of g = 0, £ = 100, we find y/ (x 2 )/ ~ 10 16 GeV, which is comparable to the case 
of j ~ 3 x 10~ 4 ,£ = 0. The unique feature of ^-resonance is the existence of £77 suppression effect. As £ increases 
up to about 100, (% 2 ) / also increases because %-field stays in the broad resonance bands longer. However, when £ 
exceeds about 100, £77 suppression effect by the production of x-particles is significant and (x 2 ) / docs not increase 
for the case of £ > 100. Rather, beyond £ ~ 200, (x 2 )/ decreases by £77 suppression effect as (x 2 )/ oc £~ 3 / 2 . In the 
case of g = 0, £ > 0, we find that the maximal value of (x 2 ) / is \[W) max ~ 10 16 GeV at £ = 100 - 200. 

As for the case of g = 0,£ < 0, the relation between Ak and q becomes Ak = —2q/3 + ft 2 /a 2 and the resonance 
band is further broader than the case of g = 0, £ > 0. The important difference from other cases is that A^-q curve 
will pass through the zeroth instability band below the curve of Ak = —q 2 /2- As a result, even if q decreases under 
unity by expansion of the Universe, the modes close to k — always stay in the resonance band. In this case, we 
find the termination in the growth of (x 2 ) at which the growing rate \i of (Y 2 ) balances the expansion rate of the 
Universe. In — 1 < £ < case, the increase of (x 2 ) is not expected because the Hubble expansion rate surpasses the 
growing rate. In — 3 < £ < — 1 case, although (x 2 ) increases as the passage of time, it takes much time to reach 
its maximum because fi is very small and resonance terminates before £77 suppression effect becomes significant. For 
£ < — 3, (x 2 ) / takes rather large values and (x 2 ) / takes maximal value \J (x 2 ) max ~ 2 x 10 17 GeV when £ w —4. As 
£ decreases from £ < — 4, although the growth rate \x increases, the final fluctuation of x-particle decreases. This is 
because £77 suppression effect due to x-particle production plays a crucial role to terminate the resonance. As a result, 
(x 2 )/ decreases as |£|~ 3//2 for |£| >• 1. We should also note that the value of \J (x 2 ) max m tne case of g = 0,£ < is 
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greater than that in the case of g = 0, £ > 0. This is because the resonance bands for £ < are broader than those 
for £ > 0. 

We have also studied the combined resonance of interactions g 2 4> 2 x 2 an d £,Rx 2 ■ The structure of resonance is quite 
different depending on two parameters of g and £. What we had been interested in is whether ^-coupling assists 
^-resonance in any cases. However, we find that this is not the case. For g > 3 x 10~ 4 , and back reaction effect on 
the inflaton field and metric is important, ^-resonance is not assisted by ^-coupling because of £77 suppression effect as 
well as back reaction effect at the final stage of preheating. In the parameter range where x-particles are significantly 
produced only by g-resonance, the maximal value of (x 2 ) / does not increase even if we include the ^-coupling. On the 
other hand, in the case of g < 3 x f ^-coupling may assist g-resonance in the parameter ranges of 100 < £ < 200 
and —00 < £/g 2 < — (Mp^/m) 2 /16tt. In particular, for g < 1 x 10~ 5 , the structure of resonance is almost the same 
as the case with g = 0, £ 7^ 0, and it is essentially the ^-resonance. We find that the maximal value of y/ (x 2 )/ is 
about 2 x 10 17 GeV for g < 1 x 10~ 5 ,£ w —4, which is larger than the minimally coupled case with jwlx 10~ 3 . 

There are several things we did not investigate in this paper. One of them is the rescattering effect. As ^-particles 
are produced significantly, the fluctuations of the inflaton field are also generated and would affect the production 
of x-particles. Although it is expected that the structure of resonance does not change so much at the first stage 
of preheating, the rescattering effect will modify the final production of x-particles because it becomes important 
at the final stage of preheating when x-partides are significantly produced. As (/>-particles are created, spatially 
inhomogeneity of the </>-field would prevent from the resonant production of x-particles, resulting r/ — £k 2 (x 2 ) to be 
reduced. Then we may find either insufficient production of x-particles or the same amount of x-particles with a delay 
of time to reach the limiting value due to the £77 suppression. As for </>-particle creation through the present coupling, 
we did not investigate details at present. For complete study of preheating, however, we should consider the growth 
of ^-particles due to rescattering effect quantitatively beyond mean field approximation. 

We did not consider the case when the coupling g 2 between <p and x is negative, although it was pointed out in p{}] 
that the production of ^-particles increases significantly compared with the positive coupling case. This is similar to 
the geometric reheating with negative £ in the sense that A/~ can take a negative value in the Mathieu chart. It may 
be interesting to study the combined resonance of negative g 2 and £. 

In this paper, we have also not studied the metric perturbation in preheating phase. However, several authors 
pointed out that metric perturbation is influenced by the parametric resonance pOpp^]. It was recognized the Bardeen 
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parameter is well conserved quantity in reheating phase except the short period when <j> is close to zero pc|j33[] . On 
the other hand, Bassett et al. [[?2| recently found that the rapid growth of metric perturbation by negative coupling 
instability is expected and this stimulates the growth of scalar field. It is worth investigating whether the growth of 
metric perturbation enhances the fluctuation of \ field non-minimally coupled to spacetime curvature R with £ < 0. 

We have studied a parametric resonance by ^-coupling (£i?x 2 ) as wen as ^-interaction (g 2 <j> 2 x 2 )- Although we 
expect that any scalar field will couple to spacetime curvature R through quantum effects, the value of £ considered 
here may be too large. However, in other theories of gravity such as the Brans-Dicke theory |3, the induced gravity 
p5| , and the higher-curvature theories we may have different types of coupling to the spacetime curvature, which 
might give a natural mechanism for an effective resonance. These issues are under investigation. 
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TABLE I. The final values (X 2 ) f obtained by the analytical estimation and by the numerical calculation in the case of 
g — 0,£ > 0. The analytical estimation gives a good approximation to the numerical results, (X 2 ) numerical takes the maximal 
value at £ ~ 100. We also show the time tf when the resonance ceases and £rj which indicates a suppression effect. For large 
£, the suppression effect (q ~ (1 + 6^77) _2 ) becomes crucial. 



t 


*/ 


(X ) analytic 


(X } numerical 


£,T]numerical 


20 


1.85 




8.610 


8.7 x 10~ 8 


50 


2.85 




3.090 x 10 3 


1.9 x 10~ 4 


70 


2.85 




3.750 x 10 s 


4.6 x 10~ 2 


100 


2.40 


4.138 x 10 5 


9.550 x 10 5 


2.4 x 10 _1 


500 


1.31 


1.513 x 10 s 


2.109 x 10 5 


1.3 


1000 


1.31 


5.626 x 10 4 


8.318 x 10 4 


2.1 


10000 


1.31 


1.922 x 10 3 


2.630 x 10 3 


6.6 



TABLE II. The final value (x 2 )/ obtained by the analytical estimation and by the numerical calculation in the case of 
g — 0, £ < 0. For |£j > 4, both values show the same tendency that {X ) / gets large for decreasing £|. However, for |£| < 3, the 
analytic estimation by Eq. (4.8) can not be applied (see text for the detail), (x 2 ) numerical takes the maximal value at £ w —4. 





*/ 


{X ) analytic 


{X ) numerical 


|4 ^numerical 


-2 


7.5 x 10 b 


1.049 x 10 a 


3.184 x 10 3 


3.2 x 10"' 


-4 


3.3 x 10 5 


5.425 x 10 8 


3.162 x 10 8 


1.3 x 10 _1 


-10 


1.0 x 10 3 


1.758 x 10* 


1.122 x 10 8 


2.8 x 10 _1 


-20 


28.28 


6.903 x 10 7 


4.130 x 10 7 


4.2 x 10 _1 


-50 


15.26 


1.896 x 10 7 


1.099 x 10 7 


6.9 x lO" 1 


-100 


7.66 


6.994 x 10 6 


4.355 x 10 6 


1.1 


-1000 


4.61 


2.355 x 10 s 


1.606 x 10 5 


4.0 


-10000 


4.18 


7.591 x 10 3 


5.069 x 10 3 


12.7 



TABLE III. The numerical result for the broad resonance in the g = 3 x 10~ 4 case, ib and (x 2 )f> are the values when the 
broad resonance regime ends. (x 2 }fc gets larger for large £ because the resonance band becomes very broad. However, (x 2 )/ 
turns out to be smaller for large £ because of the suppression effect, although it does not depend on £ so much. 
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<x 2 > 6 


if 


<x 2 >/ 


-20 


6.26 


3.019 


14.36 


2.951 x 10 b 





6.26 


5.495 


14.36 


7.244 x 10 5 


50 


6.26 


3.162 x 10 


12.85 


6.309 x 10 5 


100 


6.26 


2.344 x 10 2 


11.88 


4.570 x 10 s 
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TABLE IV. The numerical result for the new-broad resonance on the A k — 2q/3 line on the Mathieu chart, g does not take 
large value for this parameter range of £. The effective resonance is expected only for |£| > 200. 



t 


9 


(X 2 )f 


-20 


2.745 x 10" b 




-50 


4.341 x 10~ 5 




-70 


5.137 x 10~ 5 




-100 


6.140 x 10~ 5 


5.012 x 10 


-200 


7.089 x 10~ 5 


1.013 x 10 5 



TABLE V. The numerical result for the wide resonance on the Ak = line on the Mathieu chart. The resonance occurs 
mainly by the ^-coupling. The effective resonance is expected only for |£| > 100. 





9 


(x 2 )/ 


-20 


2.241 x 10" 5 




-50 


3.545 x 10~ 5 


3.981 


-70 


4.194 x 10~ 5 


1.016 x 10 4 


-100 


5.013 x 10~ 5 


3.206 x 10 5 


-200 


7.089 x 10~ 5 


1.047 x 10 6 
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Figure Captions 



FIG. 1: The schematic diagram of the Mathieu chart and the typical paths for three 

types of resonance. The lined regions denote the instability bands (zeroth, first, 
second, - • •). The line of Ak = 2q is the typical line of the ordinary (/-resonance, 
while the lines of Ak = 2q/3 and Ak = — 2g/3 show the lowest limits of the 
resonances by positive and negative ^-couplings, respectively. We find that the 
width of the instability bands becomes wider for large q. The Floquet index in 
the lower instability band gets larger for fixed q. 

FIG. 2: The evolution of (X 2 ) as a function of i in the case of g = 0, £ > ((a) £ = 50, 

(b) £ = 100 and (c) £ = 1000). We find that the parametric resonance occurs 
by the positive ^-coupling, (X 2 ) increases exponentially and reaches to its final 
value (x 2 )/- After then, it decreases by the adiabatic expansion, (X 2 )/ takes 
the maximal value for £ ~ 100. For larger £ (£ > 100), although the growth 

rate becomes large, the final value (X 2 )/ is suppressed by the £77 term. 

FIG. 3: The final value of (X 2 ) as a function of £ in the case of g = 0, £ > 0. For 

^ £ ^ 1^0' (x 2 )/ increases as £ increases because the initial value of q gets 
large. For 100 < £ < 200, (X 2 )/ takes its maximal value (~ 10 6 ), and for 
£ > 200, (X 2 )/ decreases by the suppression factor as <~ £~ 3 / 2 . 

FIG. 4: The evolution of (X 2 ) as a function of t in the case of g = 0, £ = 200, m x = m. 

The significant production of massive x _ P ar ticles of order m ~ 10 13 GeV is 
expected by the positive ^-resonance. 

FIG. 5: The evolution of (X 2 ) as a function of t in the case of g — 0, £ < ( (a) £ = — 20 

, (b) £ = —50 and (c) £ = —100). (X 2 ) increases exponentially in the first stage 
by the parametric resonance and reaches to its final value (X 2 )/- After then, 
it approaches constant value because the creation rate of x-particle and the 
expansion rate of the Universe balance. 

FIG. 6: The evolution of (X 2 ) as a function of i in the case of g — 0, £ = —200, m x = m. 

The final value of (X 2 ) / ~ 10 6 5 is larger than the case of g — 0, £ = 200, m x = 
m. 

FIG. 7: The final value of (X 2 ) as a function of fh x = m x /m in the case of £ = —1000. 

(X 2 ) / decreases as m x increases, and parametric resonance can not be expected 
for fh x > 100 (namely, for m x > 10 15 GeV). 

FIG. 8: The relation between Ak and q in the case of the combined resonance of g and £. 

(a) As £/g 2 decreases from 00 to — (MpL/m) 2 /247r, the relation changes from 
Ak = 2g/3 + k 2 /a 2 to q = as the arrow in the figure. We call the resonance 
in each parameter region as follows. < £/g 2 < 00 : new- broad resonance, 
and — (Mpl/to) 2 /247t < £/g 2 < : ordinary-broad resonance, (b) As £/g 2 
decreases from — (Mpl/to) 2 /247t to —00, the relation changes from q = to 
Ak = —2/'5q + k 2 /a 2 as the arrow in the figure. We call the resonance in each 
parameter region as follows. — (Mpl/to) 2 /167t < £/g 2 < — {Mp^ / m) 2 / 2Att 
: ordinary-broad resonance, — (MpL/m) 2 /127r < £/.g 2 < — (Mp^ / m) 2 / 16ir : 
new-broad resonance, and —00 < £/<? 2 < — (MpL/m) 2 /127r : wide resonance. 

FIG. 9: The evolution of (X 2 ) as a function of t in the case of g = 3.0 x 10~ 4 ( (a) £ = 0, 

(b) £ — 50 and (c) £ = 100). We find that the resonance is divided into two 
stages; one of which is the first broad resonance stage and the other is the 
narrow one. The growth rate in the broad resonance stage becomes larger for 
larger £. However the final value (X 2 )/ is suppressed by the £77 term. 
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FIG. 10: The evolution ol (x 2 ) as a function of i in the case of g = 1.0 x 10~ 3 ( (a) £ = 0, 
(b) £ = 50 and (c) £ = 100). The broad and narrow resonance stages are not 
distinguished. The growth rates are almost the same for each £, because the 
resonance occurs mainly by (/-coupling. When (x 2 ) reaches its maximal value, 
the back reaction effect terminates the resonance. 

FIG. 11: The structure of resonance, the back reaction and £77 suppression effects in 
terms of £ and g. The regions [A], [B], [C], [D] denote new-broad resonance 
(£ > 0), ordinary-broad resonance (— (Mp^ / m) 2 / 16ir < £/g 2 < 0), new-broad 
resonance (— (Mp^/m) 2 /12ir < £/g 2 < — (Mpi, / m) 2 / 16n) , wide resonance 
(—00 < £/g 2 < — (Mpl / m) 2 / 12tt) respectively. With this diagram, we eas- 
ily understand the resonance structure. The lined regions (g > 3 x 10~ 4 and 
g > 5 x 10~ 6 |£|) denote those where the back reaction effect is significant. The 
shaded regions (g < 3 x 10~ 4 and cither £ > 100 or £/g 2 < - 0.1(M pl /to) 2 ) 
denote those where £77 suppression effect becomes important. 

FIG. 12: (x 2 )/ m terms of £ and g. We find three plateaus, one of which (x 2 )/ takes 
maximal value {x )max ~ 3.2 x 10 8 for g < 1 x 10 _5 ,£ w -4. 
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